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ABSTRACT 


Thermoelectric  Cooler  Design  is  a  Microsoft  Windows  program  to  be  used  as  an  aid  in 
the  design  of  thermoelectric  cooler  devices.  This  program  was  written  to  be  used  to 
quickly  model  and  compare  alternative  designs.  A  couple's  optimum  coefficient  of 
performance  and  maximum  heat  pumping  can  be  quickly  determined.  Other  major 
features  of  the  program  include  the  ability  to  change  material  properties  and  dimensions 
of  couples,  analyze  cascaded  couples,  and  graph  performance  parameters.  A  brief 
description  of  thermoelectric  cooler  theory,  modelling  assumptions  and  complete  source 
code  listing  is  included. 
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I.  INTRODUCTION 


A.  FEATURES  OF  THERMOELECTRIC  COOLERS 


The  unique  features  of  thermoelectric  coolers  include: 

1 .  Compact. 

2.  No  moving  parts. 

3.  The  coefficient  of  performance  is  independent  of  capacity. 

4.  The  cooling  rate  is  easily  controlled. 

5.  Reversible  operation  is  possible  allowing  heating  as  well  as  cooling. 

6.  Vibration-free  operation. 

7.  Precise  temperature  control. 

B.  SPACEBORNE  APPLICATIONS 

Thermoelectric   coolers   offer   an   alternative   solution   to   cooling   and   heating 
requirements  of  spacecraft.  Spaceborne  applications  for  thermoelectric  coolers  include: 

1.  Equipment  cooling  —  removal  of  waste  heat  in  order  to  maintain 
components  at  their  optimal  operating  temperature.  Thermoelectric  coolers 
are  especially  useful  for  the  control  of  small  hot  spots. 


2.  Life  support  temperature  control  —  thermoelectric  coolers  are  very  reliable 
making  them  ideally  suited  for  long  duration  missions  such  as  the  space 
station  or  manned  Mars  mission. 

3.  Cooling  of  infrared  detectors  —  thermoelectric  coolers  eliminate  the  need 
for  coolant  and  associated  support  equipment.  Of  even  more  importance  is 
the  elimination  of  consumable  coolants  which  also  increases  the  life  of  the 
spacecraft. 

4.  Microgravity  applications  —  thermoelectric  coolers  offer  precise, 
vibration-free  temperature  control  for  free-flying  orbiting  microgravity 
research  platforms. 

Reliability,  the  lack  of  moving  parts  or  consumables,  and  vibration-free  operation 
are  especially  important  to  spaceborne  cooling  systems. 


II.  THERMOELECTRIC  COOLER  THEORY 


A.  THERMOELECTRIC  COUPLE  CONFIGURATION 

In  it's  most  basic  form,  a  thermoelectric  couple  consists  of  two  dissimilar  elements 
with  contact  at  two  junctions.  Figure  1  shows  a  simple  thermoelectric  circuit  consisting 
of  two  elements,  physically  joined  on  one  side  while  electrically  joined  at  the  other. 


Figure  1  Thermoelectric  Cooler 


B.  THERMOELECTRIC  EFFECTS 

The  background  for  this  section  was  gathered  from  Kraus  [Ref.  12:pp.  294-303]  and 
Ioffe  [Ref.  2:pp.  96-128].  Three  thermoelectric  effects  are  known;  the  Seebeck,  Peltier, 
and  Thomson  effects. 

1.  Seebeck  Effect 

Consider  the  junction  of  two  materials,  A  and  B.  If  a  temperature  differential 
exists  between  the  two  junctions,  one  junction  may  be  referred  to  as  the  hot  side  and  the 
other  the  cold  side.  An  open  circuit  voltage  can  be  measured  between  the  two  junctions. 
This  net  conversion  of  thermal  energy  into  electrical  energy  under  zero  current 
conditions  is  known  as  the  Seebeck  effect.  The  voltage  observed,  referred  to  as  the 
Seebeck  voltage,  is  proportional  to  the  temperature  differential  between  the  hot  and  cold 

sides  of  the  couple: 

dEs  «*  dT 

dEs  =  ±  adT 
Es  =  ±    J  adT   (V)  (1) 


where 


Es  =  Seebeck  voltage  (V) 

a  =  Seebeck  coefficient  of  material  (V/°C) 

In  the  case  of  the  couple  composed  of  materials  A  and  B, 


a  =    a^    +  I  OLfi 


where 

a^  =  Seebeck  coefficient  of  material  A 
QLg  =  Seebeck  coefficient  of  material  B 

2.  Peltier  Effect 

The  Peltier  effect  is  localized  at  the  junctions  between  materials  A  and  B.  When 

a  current  flows  through  the  junction,  heat  is  either  generated  or  absorbed  depending  on 

the  direction  of  the  current  flow.  The  amount  of  heat  is  proportional  to  the  current  and  is 

known  as  the  Peltier  heat: 

dQp  ~  Idt 

dQp  =  tnldt 

qp  =  ±nl    (W)  (2) 


where 


Qp  =  Peltier  heat  (J) 

qp  =  Peltier  heat  flow  (W) 

/  =  Current  (A) 

k  =  Peltier  coefficient  of  material  (V) 

3.  Thomson  Effect 

This  effect  concerns  the  absorption  or  generation  of  heat  in  a  conductor  carrying 

a  current  in  the  presence  of  a  temperature  gradient.  The  amount  of  heat  is  proportional  to 

the  current  and  temperature  gradient: 

dQt  -  IdTdt 

dQt  =  ±  al  dTdt 


dqt  =  taldT 
q,  =  ±al  \  dT     (W)  (3) 


where 


Ql  =  Thomson  heat  (J) 

qi  =  Thomson  heat  (W) 

/  =  Current  (A) 

a  =  Thomson  coefficient  of  material  (V/°C) 


Because  the  product  of  o  and  AT  is  a  voltage,  the  Thomson  voltage  is  therefore: 


T2 


Et=  ±    j  GdT   (V)  (4) 


C.  IRREVERSIBLE  EFFECTS 

1.  Joule  Effect 

The  evolution  of  heat  when  a  current  passes  through  a  material  may  be  referred 

to  as  the  Joule  heat  flow: 

qj  =  I2R    (W)  (5) 


In  the  case  of  a  couple,  the  elements  A  and  B  are  connected  in  electrical  series 
regardless  of  the  direction  of  the  current  flow: 

R  =  2d^  +  ^    (Q)  (6) 

where 


^5 


Py4' P5  =  electI*ical  resistivity  of  materials  >4  and  5    (Q.  cm) 

L^,Lfi  =  Length  of  materials^  and  B    (cm) 

A^,Ag  =  Cross  sectional  area  of  materials  A  and  B    (cm^) 

2.  Fourier  Effect 

The  conduction  of  heat  through  a  material  is  given  by: 


qf=  f^-W    (W)  (7) 


and,  in  the  case  of  a  couple  with  the  elements  in  thermal  parallel, 


where 


K  =  -2— I-  +  -S—S-     (W/°C)  (8) 

LA  LB 


k^,kg  =  Thermal  conductivity  of  materials  A  and  B    (W/cm°C) 

L^,Lg  =  Length  of  materials  ,4  and  B    (cm) 

A^,  Ag  =  Cross  sectional  area  of  materials  ,4  and  B    (cra^  ) 


D.  NET  HEAT  ABSORBED  AT  COLD  JUNCTION  OF  COUPLE 


In  this  analysis,  the  Thomson  heat  will  be  neglected  which  greatly  simplifies  couple 

design  procedure  without  any  significant  error  Kraus  [Ref.  l:p.  300]  and  Cadoff  [Ref. 

3:p.  8].  To  determine  the  net  heat  pumped  by  the  couple,  the  temperature  distribution 

along  the  elements  of  the  couple  must  be  determined.  Consider  a  single  thermoelectric 
element  as  indicated  in  Figure  2,  where  Ax  is  the  cross  sectional  area  of  the  element 
normal  to  the  flow  of  current  and  heat. 


Figure  2  Thermoelectric  Element 


Here,  the  boundary  conditions  are: 

T  =  To   at  x  =  0 
T  =  Te  at  x  =  L 


Assuming  steady  state  temperature  conditions,  there  are  three  heat  quantities  which 
pertain  to  an  energy  balance  over  the  incremental  length,  Ax: 
1 .  The  Fourier  heat  entering  at  point  x: 


qx  =  -kAx 


dT 

.  dx  _ 


2.  The  Joulean  heat  generated  in  the  increment  Ax 


72pAx 


V  =        AX 


3.  The  Fourier  heat  leaving  at  point  x  +  Ax  : 


«**--**[£] 


x+Ax 


The  energy  balance  in  the  steady  state  requires  that  the  heat  leaving  must  be  equal  to 
the  sum  of  the  heat  entering  and  the  heat  generated  in  the  increment  Ax  : 

qx  +  qj  =  qx+Ax 


Substituting  for  qx,  ^^+Ax  anc*  9j  gives 


-kAx 
or  with  some  algebra, 


dT_ 
_  dx 


/2  p  Ax 
+  —j =  -kAx 

Ax 


dT" 

L  <&  J  x+Ax 


kAX(im 


Idx  J  x+Ax  I  dx  . 


dT 


I  2  pAx 

+    A =      ° 

Ax 


Dividing  through  by  Ax  yields 


kA 


fill 

L  dx  ]X+Ax 


dT 
dx  _ 


„     A 


Ac 


'2p    „ 


And  taking  the  limit  as  Ax  — >  0  gives 


limit 
x->0 


kAx 

'dT' 
_  dx . 

x+Ax 

"dT' 

.  dx . 

X 

♦1p 

Ax 

<. 

Ax 

J 

which  can  be  adjusted  to 


kA 


d2T       I1  P 


dx1        A 


+  -Vs-  =  0 


(9) 


The  general  solution  to  this  simple  differential  equation  is  obtained  through  the  use  of 
a  double  integration: 

d2T       I2P* 


dx2 


+ 


kA 


2    +C1=0 


10 


and 


/2DX2 

—  +  Cix  +  C2  =  0 


2kA2x 


l 


Applying  the  boundary  conditions  for  the  hot  side  where  x  =  0  and  T  -  To 


and 


C?  =   -  To 


;2Dv2 
T  -  T0  +  — ~-  +  Ca  =  0 


2k  A- 


1 


At  the  cold  side  of  the  element  where  x  =  L  and  T  =  Te'. 


2       2 
Te  -  To  +  — £__.  +  CjZ,  =  0 

2^v 


The  solution  for  C  j ,  the  first  constant  of  integration,  is  seen  to  be 


Cl=        L 


To-Te       IzpL 


2k  A- 


so  that  the  particular  solution  of  this  simple,  second  order,  differential  equation  which 
gives  the  temperature  profile  as  a  function  of  distance  along  the  element: 


ikAi 


f  t2     t^ 
\UA2X  j 


x  =  0 


or 


11 


T(x)  =  - 


f  t2 

(£P_ 

2kA2XJ 


( 


x2- 


To  -Te 


2kA2XJ 


x  +  To 


(10) 


which  is  seen  to  be  parabolic. 


Figure  3  illustrates  a  temperature  distribution  for  a  typical  element  with  the  following 
parameters: 


z  =  3.08xl0"3     (°C    l) 
R  =  1.73xl0-3    (Q) 
K  =  3.64xl0"2    (W/°C) 
a  =  4.4X10-4    (V/°C) 
Tc  =  308.6    K 
AT  =  50    °C 


(11) 
(12) 
(13) 
(14) 
(15) 
(16) 


Due  to  the  parabolic  temperature  distribution,  the  Joulean  heat  will  be  transferred  to  both 
the  cold  and  hot  junctions  of  the  couple. 

Differentiating  equation  (10),  setting  the  derivative  equal  to  zero, 


dT_  I2pX       To  -Te       I2?  L 

dx   '  K,2  L  ^_*2 


kA 


2k  A 


permits  a  solution  for  x  =  X  which  is  the  position  where  the  maximum  temperature 
occurs: 


IzpL      To-Te 


X 


2k  A  ■ 


I2pX 

kA2 


12 


c 
£ 

0)  O 

°  O 

<      £ 

o  m 

5  »- 

5  2 

a. 
£ 

H 

ro 

(U 

I 


X  ajniejsdaiaj. 
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Rearranging  and  substituting  AT  for  To  -  Te'. 

x=  L       kATAJ 

2  I2pL 

Assuming  x  =  L  and  defining  /  as  the  fraction  of  the  Joulean  heat  transferred  to  the 
cold  junction: 

f=  L^X     1,  kATA* 

J  L  2        /2pI2 

or 

where  K  and  R  are  defined  by  equations  (6)  and  (8),  respectively. 

The  net  heat  absorbed  at  the  cold  junction  of  the  couple  is  equal  to  the  summation  of 
the  Peltier  heat,  minus  the  fraction  of  Joulean  heat  transferred  to  the  cold  junction: 


qnet  =  aITc-fI2R 


and  upon  substitution  of  the  expression  for  /  given  by  equation  (17): 


qnet  =  aITc  -\l2R-  KAT   (W)  (18) 


Therefore  the  net  heat  pumped  is  the  Peltier  heat  developed  at  the  cold  junction  less 
the  sum  of  one  half  the  Joulean  heat  produced  in  the  couples  and  the  heat  due  to 
conduction  from  the  hot  junction  to  the  cold  junction. 
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III.  THERMOELECTRIC  COOLER  MODEL 


A.  INTRODUCTION 

The  equations  in  this  chapter  were  used  to  model  the  thermoelectric  cooler  couple 
operating  parameters  and  were  given  by  Kraus  [Ref.  l:pp.  303  -  326]. 


B.  HEAT  PUMPING 

1.  Current  Yielding  Maximum  Heat  Pumped 

The  equation  for  net  heat  pumped 


qnet  =  alTc  ~  \l2R~  KAT   (W)  (18) 


may  be  differentiated  with  respect  to  current,  /.  When  the  derivative  is  set  equal  to  zero, 
the  current  that  yields  maximum  heat  pumped  is  obtained: 


%  =  aTc  -IR  =  0 

al 


15 


Solving  for  /: 


lm  =  ^p-    (A)  (19) 


2.  Maximum  Heat  Pumped 

Substituting  lm  into  the  equation  for  the  net  heat  pumped  yields  the  maximum 
heat  pumped: 

a2  T2 
qm  =  -j^~  -  KAT   (W)  (20) 


C.  MAXIMUM  TEMPERATURE  DIFFERENTIAL 

Heat  pumped  as  a  function  of  current  is  shown  in  Figure  4  using  the  parameters 

defined  in  equations  (11)  through  (16).  Setting  equation  (20)  which  is  the  maximum  heat 
pumped  equal  to  zero  and  solving  for  AT,  one  obtains  the  maximum  temperature 
differential  between  the  hot  and  cold  sides  of  the  couple: 


a2T2 
ATm  =  ~2KR      (K)  (21) 


D.  OPTIMUM  RATIO  OF  ELEMENT  CROSS  SECTIONAL  AREAS 

The  maximum  temperature  differential  given  by  equation  (21),  can  be  increased 
by  minimizing  the  product  KR  from  equations  (6)  and  (8): 


16 
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KR  =  [VALA   +  VBLB\  (kAAA   +  Mj> 


(22) 


and  if  LA  =  LB   an  expansion  gives 


KR  =  kApA  +kBpA-?-  +  kApBjt-+kBpB  (23) 


An  optimization  of  equation  (23)  can  be  accomplished  by  finding  the  point 
where  the  derivative  with  respect  to  the  area  ratio,  AA IAB ,  vanishes: 


Jiff^  =  *APB  -  hPA{%)2  =  ° 


(24) 


Then,  a  solution  of  equation  (24)  yields  the  area  ratio  that  minimizes  KR.  Thus,  in  turn 

yields  the  optimal  ATm  '■ 

A  a         \pA kR 

(25) 


AB        \PBkA 

The  area  ratio  given  by  equation  (25)  when  substituted  into  equation  (22),  yields 
the  optimum  value  of  KR,  designated  by  <j) 


or 


$  =  kAPA+2  JkAPAkBPB   +  kB ?B 


<t>  ■  (jkA?A    +  J^Pb)  <26> 
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Equation  (16)  has  the  restrictions  that  L^  -  Lg  and  that  the  area  ratio  be  determined 
using  equation  (25). 


E.  FIGURE  OF  MERIT  OF  ELEMENT  MATERIALS 

Substituting  <j>  ,  equation  (26), into  the  equation  for  Aim   equation  (21) 


a2  ^2 


A7m  =  ^tTc 


2<t> 
or 


ATm  =  —  zTc 


1_     rrl 

2 


where  z,  the  figure  of  merit  of  the  materials,  is  defined  as: 


..  a2 

which  may  be  rewritten  as: 

oi 

Z  =   


>2 
JkAPA  +JkBPB 


When  L^  =  Lg  and  the  area  ratio  determined  from  equation  (25),  the  figure  of 
merit  is  a  function  of  the  properties  of  the  elements.  To  maximize  the  figure  of  merit,  the 

materials  used  should  have  a  high  Seebeck  coefficient,  low  thermal  conductivity,  and  low 

electrical  resistivity  Cadoff  [Ref.  3:  p.  21]. 
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F.  POWER  REQUIRED  TO  OPERATE  THE  COUPLE 

To  operate  the  couple,  a  voltage  equal  to  the  sum  of  the  Seebeck  voltages  in  the 
couple  and  the  resistive  voltage  drop  must  be  applied  Kraus  [Ref.  1:  p.  304].  The  power 
required  to  operate  the  couple  is  then: 

p  =  IV  =  /(a  AT  +  IR)     (W)  (27) 


G.  COEFFICIENT  OF  PERFORMANCE 

The  coefficient  of  performance,^  ,  measures  the  performance  of  the  couple  and  is  the 
heat  pumped,  q,  divided  by  the  power  required  to  operate  the  couple,  P: 


1       P  (28) 


Substituting  equation  (18)  for  q  and  equation  (27)  for  P  gives 

aITc  -h2R-KAT 
^  =         I(IR  +  aA7)  (29) 

1.  Coefficient  Of  Performance  At  Maximum  Heat  Pumping 

One  may  express  qm  in  terms  of  AT  and  ATm  using  equations  (20)  and  (21): 
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a2^  £a2rg       „Ar 


?w  =  £Arm  -  a:aj  =  ^(Arm  -  ad 


i  -  ^-1 


In  similar  fashion,  Pm   may  be  expressed  in  terms  of  the  temperatures  by 
substituting  equation  (19)  into  equation  (27): 

Pm  =  aImAT+  ImR 


a2TcAT      a2T% 


R  R 


<*2T2C  (AT      . 
R      {Tc 


IKo^tI  fl  +  L 


2KR      V         Tc 


=  2KATm  I  I  +  4r  I  (31) 


Tc 


To  obtain  the  coefficient  of  performance  when  the  couple  is  adjusted  for 
maximum  heat  pumping,  substitute  equations  (30)  and  (31)  into  equation  (28)  for  the 
coefficient  of  performance: 
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qm        KATm(\-AT/ATm) 
^m  "  Pm        2KATm(\+AT/ATc) 

\-ATIATm 


Am  = 


2(\+AT/ATc) 


2.  Current  Yielding  Optimum  Coefficient  Of  Performance 

The  expression  for  the  coefficient  of  performance,  equation  (29), 

aITc  -h2R-KAT 
^  =         1(1  R  +  a  AT) 

may  be  differentiated  to  find  the  current  that  yields  the  maximum  coefficient  of 
performance.  This  current  is  the  value  of/  that  causes  dx\l  dl  to  vanish: 


dr\  _  I(IR  +  aAT)(aTc-IR)-(aITc-I2R/2- KAT) (2IR  +  aAT) 


or 


or 


dl  fiilR  +  aAT)1 

(I2R  +  aIAT)(aTc  -  IR)  =  (aITc  -  I2R/2  -  KAT)(2IR  +  a  AT) 

Expansion  gives: 

I2RaTc  -  I3R2  +  a2ITcAT  -  aI2RAT  =  2I2RaTc  +  a2ITcAT 
-I3R2  -  aI2RAT/2  -  2IRKAT  -  aKAT2 

aI2RAT/2  +  aTcI2R  -  2KIRAT  -  KaAT2  =  0 
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which  is  actually  a  quadratic  in  /: 


Ra{^2~  +  Tc)j2  "  ^2KRAT^1  ~  (ku&T2)  =  0 


which  possesses  a  solution 


2KRAT±  J(2KRAT)2  +  4 KR a2 AT2 [(AT/ 2)  +  Tc] 
t  =  * . LI . —  (-ir\ 

2Ra[{AT/2)  +  Tc]  V     } 


The  value  of  the  radical  in  equation  (32)  will  always  be  greater  than  2KRAT 
because  it  contains  the  sum  of  (2KRAT)2  and  another,  positive  term.  For  cooling  to 

occur,  the  current  must  not  be  negative  so  the  solution  in  which  the  radical  term  is 
subtracted  is  ignored: 


2KRAT+  J(2KRAT)2  +  4KRa2AT2[(AT/2)  +  Tc] 
''  2Ra[(AT/2)  +  Tc] 


factoring  4  (KR)2  AT2  from  under  the  radical  gives: 

2KRAT+  2KRATj\+  (a2/KR)  [(AT/2)  +  Tc] 


I  = 


2Ra[(AT/2)  +  Tc] 


or 


1  +  J\+  z[(AT/2)  +  Tc] 
1  "  ^Ar         Ra[(AT/2)  +  Tc] 
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Because 


^-  +  Tc  =  -2-= +  Tc  =  -Z-z 

2  2  2 


we  define  the  average  temperature  may  be  defined  as: 


f=   Th  +  Tc 


(33) 


and 


1         a 


1  +  JT+  z[(Th  +  Tc)/2] 


(Th  +  Tc)/2 


(34) 


Multiplying  equation  (34)  by 


Ra 
Ra 


(JT7zJ-l] 


=  1 


J\+zT -1 
and  substitution  of  the  average  temperature  given  by  equation  (33)  one  obtains: 


/  = 


KR 


OlAT]   1    +   yll+zT 


Ra2 


Ji 


+  zT  -1 


JT+zT  -1 


With  \lz  =  KR  /a   ,  expansion  gives: 


/  = 


zT 


U 


Tj\+zT -1 
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Rearranging  yields  the  expression  for  the  current  that  optimizes  the  coefficient  of 
performance: 

lo  = 


a  AT 


R  fl 


+  zT  -1 


(35) 


Figure  5  is  a  plot  of  the  coefficient  of  performance  as  a  function  of  current  based 
on  the  conditions  specified  in  equations  (11)  through  (16). 

3.  Optimum  Coefficient  Of  Performance 

Substituting  lo  into  the  coefficient  of  performance  equation  (35)  yields  the 
optimum  coefficient  of  performance: 


i\o  = 


AT 


Jl+zT-    -± 
j_c_ 


+  zT  +  1 


(36) 


As   z   approaches   <»,  x\o    approaches   the   ideal,    or   Carnot,    coefficient   of 
performance  which  is  the  term  to  the  left  of  the  brackets  in  equation  (36)  Egli  [Ref.  4:  p. 

31]  and  Tipler  [Ref.  5:pp  575-576]. 


H.  CASCADED  COOLERS 

If  a  large  AT  is  required,  a  multistage  or  cascaded  cooler  can  be  used.  In  a  cascaded 
cooler  the  heat  rejected  by  the  first  stage  of  couples  is  the  heat  load  fed  to  the  second 

stage.  This  arrangement  can  be  extended  to  any  number  of  stages.  For  best  results  when 

designing  a  multistage  cooler  for  the  lowest  temperature  attainable,  each  succeeding 
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stage  should  require  less  power  than  the  previous  stage.  It  is  also  possible  to  increase  the 
coefficient  of  performance  when  the  required  temperature  drop  is  close  to  the  maximum 
for  a  given  figure  of  merit  Ioffe  [Ref.  2:  p.  116]. 

1.  Overall  Heat  Pumped  By  Cascaded  Coolers 

The  heat  pumped  by  a  battery  of  n- stages  of  cascaded  coolers  is  equal  to  the 
summation  of  the  heat  load,  q  j  at  the  first  stage  and  the  power  required,  P,  of  each  stage 

over  all  stages  in  the  battery: 

N 

1 overall    =  a\+  2  Pf    (W) 
i  =  l 

2.  Overall  Coefficient  Of  Performance 

The  coefficient  of  performance  of  a  battery  ^-stages  of  cascaded  coolers  is  the 
heat  load  of  the  first  stage,  q  j ,  divided  by  the  total  power  required  to  operate  the  battery 

of  coolers. 

?1 
^overall  "     ^ 

I    Pi 

i=l 

I.  SUMMARY  OF  SIMPLIFYING  ASSUMPTIONS 

1.  Thomson  Voltage 

Neglect  of  the  Thomson  voltage  results  in  no  significant  error.  However,  designs 
utilizing  this  model  will  be  conservative,  understating  couple  performance  Kraus  [Ref. 
l:p.  322]. 
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2.  Temperature  Losses 

Temperature  losses  between  the  thermoelectric  couple  and  the  heat  load  or  heat 
exchanger  are  not  taken  into  account  directly.  However,  once  determined,  these 
temperature  losses  can  be  compensated  for  by  adding  the  drop  to  the  cold  side 
temperature  and  subtracting  the  drop  from  the  hot  side  temperature. 

3.  Parallel  Thermal  Paths 

Insulation  between  couples  is  assumed  to  be  perfect,  eliminating  parallel  thermal 
path  heat  losses. 

4.  Junction  Resistance 

Junction  resistance  is  assumed  to  be  zero. 

5.  Figure  Of  Merit 

In  order  to  maximize  the  figure  of  merit,  the  couple  material  lengths  are  assumed 
to  be  equal  and  area  ratio  determined  by  equation  (25). 

6.  Resistivity  And  Conductivity 

Electrical  resistivity,p  ,  and  thermal  conductivity  ,  k,  are  assumed  to  be  constant 
with  respect  to  temperature. 

7.  Couple  Power  Supply 

The  couple  power  supply  is  assumed  to  be  direct  current. 
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IV.  THERMOELECTRIC  COOLER  DESIGN  PROGRAM 


A.  INTRODUCTION 

1.  Features 

Thermoelectric  Cooler  Design  is  a  Microsoft  Windows  program  to  be  used  as  an 
aid  in  the  design  of  thermoelectric  cooler  devices.  Use  of  the  Microsoft  Windows 
graphical  user  interface  allows  the  user  who  is  familiar  with  the  Windows  environment  to 
concentrate  on  the  design  process  instead  of  a  new,  nonstandard  interface.  This  program 
was  written  to  be  used  to  quickly  model  and  compare  alternative  designs.  Modifications 
to  a  design  can  be  easily  incorporated.  Major  features  of  the  program  include: 

1 .  Ability  to  change  material  properties  and  dimensions  of  couples. 

2.  Determine  optimum  coefficient  of  performance  of  each  couple. 

3.  Determine  maximum  heat  pumping  capacity  of  each  couple. 

4.  The  ability  to  analyze  cascaded  couples. 

5.  Graph  the  performance  parameters  of  couple. 

6.  Provide  on-line  help  in  the  form  of  a  glossary,  diagrams,  and  a  temperature 
conversion  calculator. 
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2.  System  Requirements 

The  minimum  software  and  hardware  requirements  needed  by  the  computer 
system  to  run  Thermoelectric  Cooler  Design  are: 

1.  Microsoft  Windows,  version  3.0  or  later. 

2.  A  hard  drive  with  500  kilobytes  of  free  space. 

3.  A  mouse  that  is  supported  by  Windows. 

3.  Disk  Contents 

Files  included  on  the  installation  disk: 

1.  VBRUN100.DLL   -   Visual   Basic  kernel   which  must  be   in   the   same 
directory  as  THERMAL.EXE. 

2.  THERMAL.EXE  -  Thermolelectric  Cooler  Design  executable  file. 

3.  README.TXT  -  This  manual  in  ASCII  format. 

Both  VBRUN100.DLL  and  THERMAL.EXE  must  be  installed  in  order  to  run 
Thermoelectric  Cooler  Design. 

4.  Installation 

To  install  Thermoelectric  Cooler  Design: 

1 .  Insert  the  installation  disk  in  A:  drive. 

2.  Exit  to  DOS.  (Ensure  C:\  prompt  is  present). 

3.  Type  MD  THERMAL  then  press  the  ENTER  key. 

4.  Type   COPY   A:\*.*   C:\THERMAL   then   press   the   enter   key.    These 
commands  create  a  new  directory  named  THERMAL  on  the  C  drive  and 
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copy  the  files  from  the  installation  disk  into  the  C:\THERMAL  directory. 
However,  this  is  only  one  of  many  different  ways  to  copy  the  files  from  the 
installation  disk  to  the  computer.  If  the  choice  is  to  install  Thermoelectric 
Cooler  Design  on  another  drive  or  in  another  directory,  both 
VBRUN100.DLL  and  THERMAL.EXE  must  be  installed  together  in  the 
same  directory.  The  program  will  not  run  if  this  is  not  the  case. 

5.  Run  Microsoft  Windows. 

6.  Select  the  Program  Manager. 

7.  Open  the  Group  Window  in  which  Thermoelectric  Cooler  Design  is  to  be 
installed. 

8.  Choose  the  option  "new"  from  the  file  menu  of  the  Windows  Program 
Manager  window.  The  new  program  object  dialog  box  will  then  appear. 

9.  Select  the  program  item  and  choose  OK.  The  program  item  properties 
dialog  box  appears. 

10.  Type  a  description  to  appear  under  the  program  icon. 

1 1 .  Select  the  browse  button  to  display  a  list  of  files  and  directories. 

12.  Select  THERMAL.EXE  and  choose  OK. 

13.  Select  OK  to  close  the  program  item  properties  dialog  box. 

14.  The  Thermoelectric  Cooler  Design  program  is  now  installed. 

B.  USING  THE  PROGRAM 

1.  Starting  And  Exiting  The  Program 

To  run  Thermoelectric  Cooler  Design,  simply  double-click  the  program  icon. 
After  displaying  the  introduction,  the  program  menu  bar  will  appear.  This  program  menu 
bar,  titled  "Thermoelectric  Cooler  Design"  provides  access  to  all  program  functions.  To 
increase  readability,  select  the  minimize  on  use  option  under  the  Windows  program 
manager  options  menu. 
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Thermoelectric  Cooler  Design 


Info      Glossary     Qiagrams     Material  Properties      Optimize      G_raph      Calculator      Exit 


Figure  6  Program  Menu  Bar 
To  exit  the  program,  select  exit  from  the  program  menu  bar  and  choose  yes 
when  prompted.  Once  exited,  all  current  design  data  is  lost. 

2.  Setting  Material  Properties 

Prior  to  determining  the  optimum  coefficient  of  performance  or  maximum  heat 
pumping  capacity  of  a  couple,  the  material  properties  and  dimensions  of  the  couple  must 
be  input. 

To  begin,  select  material  properties  from  the  program  menu  bar.  The  material 
properties  menu  will  appear  with  default  values.  If  these  values  are  acceptable,  choose 
OK  and  proceed  to  the  optimize  menu.  If  the  default  values  are  not  acceptable  they  may 
be  changed  by: 

1.  Choosing  either  the  first  or  second  material.  The  a,  p  and  k  of  the  selected 
material  is  then  displayed.  All  values  in  blue  can  be  changed  by  the  user. 
Appropriate  values  of  a ,  p  and  k  may  be  entered  by  the  user  by  highlighting 
the  text  to  be  changed  with  the  mouse  and  inputting  the  desired  value  from 

the  keyboard. 

2.  Changes  to  the  dimensions  of  the  couple  are  input  in  a  similar  manner.  Both 
the  area  and  length  of  the  second  material  may  be  input  directly  by  the  user. 

3.  Once  the  desired  material  properties  and  couple  dimensions  have  been 
input,  the  overall  resistance,  thermal  conductivity,  and  figure  of  merit  are 
calculated  by  selecting  the  update  button.  In  order  to  proceed  to  the  optimize 
menu,  select  OK.  To  return  all  values  to  the  defaults,  choose  the  default 
button. 
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"  Properties 
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Resistance  (R) 
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1.719E  03 

a 
w/°c 

3.659E-02 
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Figure  7  Material  Properties  Window 


3.  Optimum  Coefficient  Of  Performance 


To  calculate  the  optimum  coefficient  of  performance  of  a  couple: 


1 .  Ensure  that  the  material  properties  menu  has  been  reviewed. 

2.  Select  the  coefficient  of  performance  option  of  the  optimize  menu. 

3.  The  performance  window  will  then  appear.  For  a  given  low  temperature  the 
maximum  high  temperature  is  computed.  The  default  A7"  is  the  maximum 
and  may  be  changed  by  manipulating  the  vertical  scroll  bar.  Once  the  heat 
load,  low  temperature  and  AT  are  acceptable,  the  optimum  coefficient  of 
performance  is  calculated  by  pressing  the  compute  button.  The  optimum 

coefficient  of  performance  and  associated  performance  values  are  then 
displayed.  Changes  to  the  heat  load,  low  temperature,  and  AT  can  be  made 
and  all  couple  performance  values  recomputed  by  choosing  restart. 
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Figure  8  Optimum  Coefficient  of  Performance  Window 


H 

OpUmum  Coefficient  t*f  Ptii&tmmm 

1  w 

Print 

A/couple 
Wcouple 
W/couple 

23.37 

qT 
qP 
qJ 
qF 

qo 

300.0 

W       Number  of  Couples 

W       Total  Heat  Rejected 

W       Coefficient  of  Performance 

W 

Power  Consumption 

229 

n  n^RR 

3  1733 

RH4  ft 

W 

1.33 

0.4696 

0  9843 

1.3941 

304.8 

w 

1.3102 

Figure  9  Coefficient  of  Performance  Results  Window 
4.  Maximum  Heat  Pumping 

To  calculate  the  maximum  heat  pumping  capacity  of  a  couple: 

1 .  Ensure  that  the  material  properties  menu  has  been  reviewed. 

2.  Select  maximum  the  heat  pumping  option  of  the  optimize  menu. 

3.  The  pumping  window  will   appear.   For  a  given  low  temperature  the 

maximum  high  temperature  is  computed  and  displayed.  The  default 
temperature  difference,  AT,  is  the  maximum  and  may  be  changed  by 
manipulating  the  vertical  scroll  bar.  Once  the  heat  load,  low  temperature,  and 
AT  are  acceptable,  the  maximum  heat  pumped  is  calculated  by  choosing 
compute.  The  maximum  heat  pumped  and  associated  couple  performance 
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values  are  then  displayed.  Changes  to  the  heat  load,  low  temperature  and  AT 
can  be  made  and  all  couple  performance  values  recomputed  by  choosing 

restart. 
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Figure  1 0  Maximum  Heat  Pumped  Window 
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Figure  Eleven  Heat  Pumped  Results  Window 

5.  Cascading  Stages 

Once  the  initial  couple  performance  values  have  been  computed  for  either 
optimum  coefficient  of  performance  or  maximum  heat  pumping,  cascading  of  multiple 
couple  stages  can  be  simulated  by  selecting  cascade.  The  cascade  button  is  not  visible 
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until  the  first  stage  couple  performance  has  been  calculated.   Each  additional   stage 

assumes  the  rejected  heat  of  the  previous  stage  at  the  heat  load  and  the  high  temperature 

of  the  previous  stage.  The  maximum  temperature  at  the  new  low  temperature  is  then 
calculated.  While  the  A T  of  the  previous  stage  is  assumed  as  a  default,  the  AT  for  every 

stage  of  the  cascade  can  be  changed  with  the  vertical  scroll  bar.  While  the  program  only 

displays  five  stages  at  a  time,  an  unlimited  number  of  cascade  stages  can  be  computed 

and  displayed.  Results  of  stages  beyond  five  are  shown  in  either  the  optimum  coefficient 

of  performance  or  maximum  heat  pumping  windows,  but  not  in  the  cascade  window. 
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Figure  Twelve  Cascaded  Stages  Window 

6.  Graphing  The  Coefficient  Of  Performance  And  Heat  Pumped 

Thermoelectric  Cooler  Design  is  able  to  graph  both  coefficient  of  performance 
and  heat  pumped  as  a  function  of  current.  These  graphs  are  useful  in  determining  the 
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effect  of  a  small  change  in  current  on  both  the  maximum  coefficient  of  performance  and 
the  heat  pumped.  The  effects  of  varying  the  temperature  differential  across  the  couple 
can  also  be  compared. 

To  graph  either  the  coefficient  of  performance  or  heat  pumped  as  a  function  of 

current: 

1.  Ensure  that  the  material  properties  have  been  reviewed  and  optimum 
coefficient  of  performance  or  maximum  heat  pumped  has  been  determined. 
The  graph  option  of  the  program  menu  bar  will  not  be  available  until  these 
steps  are  accomplished. 

2.  Select  the  graph  option  from  the  program  menu  bar.  The  graph  window  will 
appear. 

3.  Select  the  function  option  from  the  graph  window  menu  bar.  Choose  either 
coefficient  of  performance  versus  current  or  heat  pumped  versus  current. 

4.  Choose  draw  from  the  graph  window  menu  bar.  The  selected  function  is 
then  graphed  for  the  temperature  differential,  AT,  chosen  in  the  performance 
or  heat  pumped  window. 

5.  To  change  the  temperature  differential,  manipulate  the  AT  vertical  scroll 
bar  in  either  the  performance  or  heat  pumped  window  as  appropriate. 

To  superimpose  plots  of  representative  of  various  temperature  differentials  and 

functions,  select  superimpose  from  under  the  options  choice  of  the  graph  window  menu 

bar.  For  example,  to  compare  the  effect  of  varying  the  temperature  differential  on  a 
graph,  change  the  AT  vertical  scroll  bar  in  either  the  performance  or  heat  pumped 
window.  Subsequent  draw  commands  are  then  superimposed  on  the  current  graph.  To 

clear  the  plot,  deselect  superimpose  and  then  select  draw 


37 


f^ifarmanee 


Function     Draw     Options     Print 
Exit 


4.74  e+00 


COP 


0.0 


2.74  e+00 


1.55  e+02 


Results 
Delta  T 
Optimum 
Current 

K 
A 

130.00 

7.13  e+01 

Figure  Thirteen  Optimum  Coefficient  of  Performance  Graph 
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Figure  Fourteen  Maximum  Heat  Pumped  Graph 


38 


The  number  of  samples  used  to  produce  the  graph  may  be  adjusted  as  required. 
Select  the  sample  size  desired  under  the  samples  choice  under  the  options  selection  of  the 
graph  window  menu  bar. 

7.  Temperature  Conversion  Calculator 

To  convert  between  temperature  scales: 

1.  Select  the  calculator  option  of  the  program  menu  bar.  The  temperature 
conversion  calculator  will  appear. 

2.  Input  the  temperature  to  convert  from  by  highlighting  the  blue  numerical 
display  with  the  mouse  by  holding  down  the  left  mouse  button  while 
sweeping  the  cursor  from  left  to  right  across  the  calculator  display  text.  Once 
the  text  is  highlighted,  type  in  the  desired  temperature. 

3.  Select  the  temperature  scale  to  convert  "from". 

4.  Select  the  temperature  scale  to  convert  "to". 

5.  Choose  convert  and  read  the  temperature  in  the  calculator  display. 

6.  To  clear  the  display,  select  clear  display  option. 


Temperature 


Exit 


Figure  Fifteen  Temperature  Conversion  Calculator 
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8.  Diagrams  And  Glossary 

On-line  assistance  is  available  in  the  form  of  diagrams  and  a  glossary. 

Diagrams  of  a  single  couple  and  cascaded  couples  are  available  by  selecting  the 
diagram  option  from  the  program  menu  bar. 

A  glossary  of  terms  used  throughout  the  program  is  available.  Choose  the  glossary 
option  of  the  program  menu  bar. 
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V.  CONCLUSIONS 


The  biggest  drawback  of  a  thermoelectric  cooler  is  it's  inherent  low  coefficient  of 
performance,  but  new  materials  with  better  thermoelectric  properties  could  eventually 
make  thermoelectric  devices  competitive  with  compressor  systems.  While  unable  to 
compete  directly  with  compressor  based  systems,  thermoelectric  coolers  do  have  many 
unique  characteristics  making  them  useful  in  a  variety  of  applications.  For  spaceborne 
applications,  high  reliability,  vibration  free  operation,  and  no  dependence  on 
consumables  are  their  most  important  features. 

The  design  of  thermoelectric  coolers  is  tedious  without  the  aid  of  a  computer. 
Unfortunately,  the  utility  of  many  programs  written  for  DOS  are  handicapped  by  a  poor 
user  interface.  More  time  is  often  required  to  figure  out  how  to  use  the  program  than  to 
solve  the  initial  problem.  In  writing  thermoelectric  cooler  design,  the  Windows  graphical 
user  interface  was  used  to  provide  a  better  user  interface.  Windows  offers  the 
programmer  several  powerful  program  development  tools.  For  the  user,  Windows  offers 
a  familiar  working  environment  enabling  the  user  to  concentrate  on  the  task  at  hand. 
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APPENDIX 


Computer  code  listing  for  Thermoelectric  Cooler  Design  program. 

Thesis.txt 

//  Code  for  main  menu  bar 
Sub  ProgramInfoName_Click  () 

MsgBox  "LCDR  Clifton's  Naval  Postgraduate  School  Thesis" 
End  Sub 

Sub  DefinitionsName_Click  () 

Definitions. Show 
End  Sub 

Sub  TemperatureName_Click  () 

Screen. MousePointer  =  1 1 

TempCalc.Show 

TempCalc.  WindowState  =  0 

Screen. MousePointer  =  0 
End  Sub 

Sub  ChooseName_Click  () 

Material.  Show 
End  Sub 

Sub  ShowGraphName_Click  () 

Screen. MousePointer  =  1 1 

Graph.  Show 

Graph.WindowState  =  0 

Screen. MousePointer  =  0 
End  Sub 

Sub  ExitName_Click  () 

Screen. MousePointer  =  11 

Ending. Show 

Screen. MousePointer  =  0 
End  Sub 

Sub  Form_Load  () 

Main.Left  =  (Screen. Width  -  Main.Width)  /  2 

Main. Top  =  0 
End  Sub 
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Sub  Form_Click  () 

Main.Left  =  (Screen.Width  -  Main.Width)  /  2 
Main. Top  =  0 

End  Sub 

Sub  PerformanceName_Click  () 
End  Sub 

Sub  PumpingName_Click  () 

Pumping.  Show 
End  Sub 

Sub  Command l_Click  () 
Results.CoeffText.Text  = "" 
Results. NumberCouplesText.Text  =  "" 
Results.TotalHeatRejectedText.Text  = "'" 
Results.AperCoupleText.Text  = "" 
Results.  VperCoupleText.Text  =  "" 
Results.WperCoupleText.Text  =  "" 
Results.qTText.Text  =  "" 
Results.  qPText.Text  =  "" 
Results.qJText.Text  =  "" 
Results.qFText.Text  =  "" 
Results.qoText.Text  = "" 

End  Sub 

Sub  InfoName_Click  () 

MsgBox  "LCDR  Clifton's  Naval  Postgraduate  School  Thesis" 
End  Sub 

Sub  SinglePumpingName_Click  () 

Screen.  MousePointer  =  1 1 

ScrollUpdate  =  False 

Tc  =  Val(Pumping.LowTempText.Text) 
DeltaTMax  =  .5  *  z  *  (Tc  A  2) 
TMax  =  DeltaTMax  +  Tc 

Pumping. Tempscroll. Max  =  Tc  *  10  +  1 
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Pumping.Tempscroll.Min  =  TMax  *  10  -  1 
Pumping.Tempscroll.  Value  =  Pumping.Tempscroll.Min 

Pumping.TMinDisplay.Text  -  Format$(Tc,  "##,##0.0") 
Pumping.TMaxDisplay.Text  =  Format$(DeltaTMax  +  Tc,  "##,##0.0") 

Th  =  Pumping.Tempscroll. Value  /  10 
Pumping.HighTempDisplay.Text  =  Format$(Th,  "##,##0.0") 
DeltaT  =  Th  -  Tc 
Pumping.DeltaTDisplay.Text  =  Format$(DeltaT,  "##,##0.0") 

Pumping.  Show 
Pumping. WindowState  =  0 

Screen.  MousePointer  =  0 

ScrollUpdate  =  True 

End  Sub 

Sub  GlossaryName_Click  () 

Screen. MousePointer  =  1 1 

Definitions.  Show 

Definitions. WindowState  =  0 

Screen. MousePointer  =  0 
End  Sub 

Sub  CoolerName_Click  () 

Screen. MousePointer  =  1 1 

Cooler.Show 

Cooler.WindowState  =  0 

Screen. MousePointer  =  0 
End  Sub 

Sub  CascadeName_Click  () 

Screen. MousePointer  =  1 1 

Cascade. Show 

Cascade.  WindowState  =  0 

Screen. MousePointer  =  0 
End  Sub 

Sub  SinglePerformanceName_Click  () 

Screen. MousePointer  =11 
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ScrollUpdate  =  False 

Tc  =  Val(Performance.LowTempText.Text) 
DeltaTMax  =  .5  *  z  *  (Tc  A  2) 
TMax  =  DeltaTMax  +  Tc 

Performance. Tempscroll. Max  =  Tc  *  10+1 

Performance. Tempscroll.Min  =  TMax  *  10  -  1 

Performance. Tempscroll.  Value  =  Performance.Tempscroll.Min 

Performance.TMinDisplay.Text  =  Format$(Tc,  "##,##0.0") 
Performance.TMaxDisplay.Text  =  Format$(TMax,  "##,##0.0") 

Th  =  Performance.Tempscroll.  Value  /  10 
Performance.HighTempDisplay.Text  -  Format$(Th,  "##,##0.0") 
DeltaT  =  Th  -  Tc 
Performance.DeltaTDisplay.Text  =  Format$(DeltaT,  "##,##0.0") 

Performance.  Show 
Performance.  WindowState  =  0 

Screen. MousePointer  =  0 

ScrollUpdate  =  True 
End  Sub 
Sub  MaterialName_Click  () 

Screen. MousePointer  =  1 1 

alpha_A  =  .00023 
alpha_B  =  -  00021 
rho_A  =  .001! 
rho_B  =  .001! 
kappa_A  =  .017 
kappa_B  =  .0145 

Aa_over_Ab  =  Sqr((kappa_B  *  rho_A)  /  (kappa_A  *  rho_B)) 

Ab  =  .3848 

Aa  =  Aa_over_Ab  *  Ab 

dB  =  2  *  Sqr(Ab  /  pi) 

dA  =  2*Sqr(Aa/pi) 
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La  =  .3 18 

Lb  =  La 

alpha  =  alpha_A  +  alpha_B 

z  =  (alpha  A  2)  /  ((Sqr(kappa_A  *  rho_A)  +  Sqr(kappa_B  *  rho_B))  A  2) 

R  =  (rho_A  *  La  /  Aa)  +  (rho_B  *  Lb  /  Ab) 

K  =  ((kappa.A  *  Aa)  /  La)  +  ((kappa_B  *  Ab)  /  Lb) 

Material.  Option  1.  Value  =  True 

Material. AlphaText.Text  =  Format$(alpha_A,  "0.00000") 
Material.RhoText.Text  =  Format$(rho_A,  "0.00000") 
Material. kText.Text  =  Format$(kappa_A,  "0.00000") 

Material. FirstMaterial_a_Display.Text  =  Format$(Aa,  "##0.0000") 
Material.FirstMaterial_d_Display.Text  =  Format$(dA,  "##0.0000") 
Material. SecondMaterial_d_Display.Text  =  Format$(dB,  "##0.0000") 
Material.FirstMaterial_L_Display.Text  =  Format$(La,  "##0.0000") 
Material.RDisplay.Text  =  Format$(R,  "0.000E+00") 
Material.KDisplay.Text  =  Format$(K,  "0.000E+00") 
Material.zDisplay.Text  =  Format$(z,  "0.000E+00") 

Material.  Show 
Material. WindowState  =  0 
Screen. MousePointer  =  0 
End  Sub 

Material.txt 

//  Code  for  Material  Properties  Window 
Dim  Aa,  Ab,  La,  Lb  As  Double 


Sub  Form_Load  () 

Material.Left  =  (Screen. Width  -  Material. Width)  /  2 
Material. Top  =  (Screen. Height  -  Material. Height)  /  2 
End  Sub 

Sub  Materials_Click  () 
End  Sub 

Sub  Option  l_Click() 

If  Option  1  .Value  =  True  Then 
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AlphaText.Text  =  Format$(alpha_A,  "0.00000") 
RhoText.Text  =  Format$(rho_A,  "0.00000") 
kText.Text  =  Format$(kappa_A,  "0.00000") 

Else 

AlphaText.Text  =  Format$(alpha_B,  "0.00000") 
RhoText.Text  =  Format$(rho_B,  "0.00000") 
kText.Text  =  Format$(kappa_B,  "0.00000") 

End  If 
End  Sub 

Sub  Option2_Click  () 

If  Option  1  .Value  =  True  Then 

AlphaText.Text  =  Format$(alpha_A,  "0.00000") 
RhoText.Text  -  Format$(rho_A,  "0.00000") 
kText.Text  =  Format$(kappa_A,  "0.00000") 

Else 

AlphaText.Text  =  Format$(alpha_B,  "0.00000") 
RhoText.Text  =  Format$(rho_B,  "0.00000") 
kText.Text  =  Format$(kappa_B,  "0.00000") 

End  If 

End  Sub 

Sub  AlphaText_Change  () 
If  Option  1  .Value  =  True  Then 

alpha_A  =  Val(AlphaText.Text) 
Else 

alpha_B  =  Val(  AlphaText.Text) 
End  If 
End  Sub 

Sub  RhoText_Change  () 

If  Option  1  .Value  =  True  Then 
rho_A  =  Val(RhoText.Text) 
Else 

rho_B  =  Val(RhoText.Text) 
End  If 
End  Sub 

Sub  kText_Change  () 

If  Option  1  .Value  =  True  Then 

kappa_A  =  Val(kText.Text) 
Else 

kappa_B  =  Val(kText.Text) 
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End  If 
End  Sub 

Sub  Command l_Click  () 

Main.OptimizeName.Enabled  =  True 

Material. Hide 
End  Sub 

Sub  Command2_Click  () 
alpha_A  =  .00023 
alpha_B  = -00021 
rho_A  =  .001! 
rho_B  =  -001! 

kappa_A  =  .017 
kappa_B  =  .0145 

Aa_over_Ab  =  Sqr((kappa_B  *  rho_A)  /  (kappa_A  *  rho_B)) 

Ab  =  .3848 

Aa  =  Aa_over_Ab  *  Ab 

dB  =  2  *  Sqr(  Ab  /  pi) 

dA  =  2  *  Sqr(  Aa  /  pi) 

La  =  .318 

Lb  =  La 

alpha  =  alpha_A  +  alpha_B 

z  =  (alpha  A  2)  /  ((Sqr(kappa_A  *  rho_A)  +  Sqr(kappa_B  *  rho_B))  A  2) 

R  =  (rho_A  *  La  /  Aa)  +  (rho_B  *  Lb  /  Ab) 

K  =  ((kappa_A  *  Aa)  /  La)  +  ((kappa_B  *  Ab)  /  Lb) 

Option  1.  Value  =  True 

AlphaText.Text  =  Format$(alpha_A,  "0.00000") 
RhoText.Text  =  Format$(rho_A,  "0.00000") 
kText.Text  =  Format$(kappa_A,  "0.00000") 

FirstMaterial_a_Display.Text  =  Format$(Aa,  "##0.0000") 
SecondMaterial_a_Display.Text  =  Format$(Ab,  "##0.0000") 
FirstMaterial_d_Display.Text  =  Format$(dA,  "##0.0000") 
SecondMaterial_d_Display.Text  =  Format$(dB,  "##0.0000") 
FirstMaterial_L_Display.Text  =  Format$(La,  "##0.0000") 
SecondMaterial_L_Display.Text  =  Format$(Lb,  "##0.0000") 
RDisplay.Text  =  Format$(R,  "0.000E+00") 
KDisplay.Text  =  Format$(K,  "0.000E+00") 
zDisplay.Text  =  Format$(z,  "0.000E+00") 
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End  Sub 

Sub  Command3_Click  () 

Aa_over_Ab  =  Sqr((kappa_B  *  rho_A)  /  (kappa_A  *  rho_B)) 

Ab  =  Val(SecondMaterial_a_Display.Text) 

Aa  =  Aa_over_Ab  *  Ab 

dB  =  2  *  Sqr(Ab  /  pi) 

dA  =  2  *  Sqr(Aa  /  pi) 

Lb  =  Val(SecondMaterial_L_Display.Text) 

La  =  Lb 

alpha  =  alpha_A  +  alpha_B 

z  =  (alpha  A  2)  /  ((Sqr(kappa_A  *  rho_A)  +  Sqr(kappa_B  *  rho_B))  A  2) 

R  =  (rho_A  *  La  /  Aa)  +  (rho_B  *  Lb  /  Ab) 

K  =  ((kappa_A  *  Aa)  /  La)  +  ((kappa_B  *  Ab)  /  Lb) 

FirstMaterial_a_Display.Text  =  Format$(Aa,  "##0.0000") 
SecondMaterial_a_Display.Text  =  Format$(Ab,  "##0.0000") 
FirstMaterial_d_Display.Text  =  Format$(dA,  "##0.0000") 
SecondMaterial_d_Display.Text  =  Format$(dB,  "##0.0000") 
FirstMaterial_L_Display.Text  =  Format$(La,  "##0.0000") 
SecondMaterial_L_Display.Text  =  Format$(Lb,  "##0.0000") 
RDisplay.Text  =  Format$(R,  "0.000E+00") 
KDisplay.Text  =  Format$(K,  "0.000E+00") 
zDisplay.Text  =  Format$(z,  "0.000E+00") 

End  Sub 

Coinput.txt 

Sub  Command  l_Click  () 

Screen.  MousePointer  =  1 1 

Stage  =  0 

Performance. Hide 

History.Hide 

Results. Hide 

Screen. MousePointer  =  0 
End  Sub 

Sub  Form_Load  () 

PerformanceCmd.  Width  =  2415 

Performance. Left  =  Screen.  Width  -  Performance. Width 
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Performance. Top  =  Main. Height 
LowTemp  Update  =  True 
Cascading  =  False 

End  Sub 

Sub  Command2_Click  () 

Screen.  MousePointer  =  1 1 
Command2. Caption  =  "Wait" 
History.  WindowState  =  0 
LowTempUpdate  =  True 
Cascading  =  True 

Stage  =  Stage  +  1 

If  Stage  =  1  Then 
TotalPower  =  0 
End  If 

Tc  =  Val(LowTempText.Text) 
Th  =  Val(HighTempDisplay.Text) 

DeltaT  =  Th  -  Tc 

Th_plus_Tc  =  Th  +  Tc 

gamma  =  Sqr(l  +  (z  /  2)  *  (Th_plus_Tc)) 

Io  =  (alpha  *  DeltaT)  /  (R  *  (gamma  -  1)) 

Vo  =  (alpha  *  DeltaT)  +  Io  *  R 

P  =  Io  *  Vo 

qT  =  Val(ReferLoadText.Text) 

qP  =  alpha  *  Io  *  Tc 

qj  =  .5  *  (Io  A  2)  *  R 

qF  =  K  *  DeltaT 

qo  =  qP  -  (qJ  +  qF) 

No_of_Couples  =  qT  /  qo 

Total_Heat_rejected  =  qT  +  No_of_Couples  *  P 

Coeff_of_Performance  =  qo  /  P 

TotalPower  =  TotalPower  +  (P  *  No_of_Couples) 

If  Stage  =  1  Then 

FirstStage_qo  =  qo  *  No_of_Couples 
End  If 
OverallCOP  =  FirstStage_qo  /  TotalPower 

TMaxDisplay.Text  =  Format$(DeltaTMax  +  Tc,  "##,##0.0") 
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TMinDisplay.Text  =  Format$(Tc,  "##,##0.0") 
DeltaTDisplay.Text  =  Format$(DeltaT,  "##,##0.0") 

Results.CoeffText.Text  =  Format$(Coeff_of_Performance,  "0.0000") 
Results.NumberCouplesText.Text  =  Format$(No_of_Couples,  "#,##0.00") 
Results.TotalHeatRejectedText.Text  =  Format$(Total_Heat_rejected,  "#,##0.0") 
Results. AperCoupleText.Text  =  Format$(Io,  "#,##0.00") 
Results.VperCoupleText.Text  =  Format$(Vo,  "#,##0.0000") 
Results.WperCoupleText.Text  =  Format$(P,  "#,##0.0000") 
Results.qTText.Text  =  Format$(qT,  "#,##0.0") 
Results.qPText.Text  =  Format$(qP,  "#,##0.0000") 
Results. qJText.Text  =  Format$(qJ,  "#,##0.0000") 
Results. qFText.Text  =  Format$(qF,  "#,##0.0000") 
Results. qoText.Text  =  Format$(qo,  "#,##0.0000") 

Results.Caption  =  "Optimum  Coefficient  of  Performance"  +  "  Stage  "  +  Str$(Stage) 
Results. Show 

Select  Case  Stage 
Case  1 

History.qLoadStageOneDisplay.Text  =  Format$(qT,  "#,##0.0") 
History. qRejectedStageOneDisplay.Text  =  Format$(Total_Heat_rejected, 

"#,##0.0") 

History.  CouplesStageOneDisplay.Text  =  Format$(No_of_Couples, 

"#,##0.00") 

History. WperCoupleStageOneDisplay.Text  =  Format$(P,  "#,##0.0000") 
History.TcStageOneDisplay.Text  =  Format$(Tc,  "0.0000") 
History.ThStageOneDisplay.Text  =  Format$(Th,  "0.0000") 
History.DeltaTStageOneDisplay.Text  =  Format$(DeltaT,  "0.0000") 
History.  StageCOPStageOneDisplay. Text  =  Format$(Coeff_of_Performance, 

"0.0000") 

History.OverallCOPStageOneDisplay.Text  =  Format$(OverallCOP,  "0.0000") 
History.  Width  =  2250 
Case  2 

History.qLoadStageTwoDisplay.Text  =  Format$(qT,  "#,##0.0") 
History. qRejectedStageTwoDisplay. Text  =  Format$(Total_Heat_rejected, 

"#,##0.0") 

History. CouplesStageTwoDisplay.Text  =  Format$(No_of_Couples, 

"#,##0.00") 

History. WperCoupleStageTwoDisplay.Text  =  Format$(P,  "#,##0.0000") 
History.TcStageTwoDisplay.Text  =  Format$(Tc,  "0.0000") 
History.ThStageTwoDisplay.Text  =  Format$(Th,  "0.0000") 
History.DeltaTStageTwoDisplay.Text  =  Format$(DeltaT,  "0.0000") 
History.  StageCOPStageTwoDisplay.  Text  =  Format$(Coeff_of_Performance, 

"0.0000") 
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History.OverallCOPStageTwoDisplay.Text  =  Format$(OverallCOP,  "0.0000") 
History.  Width  =  3360 
Case  3 

History.qLoadStageThreeDisplay.Text  =  Format$(qT,  "#,##0.0") 
History. qRejectedStageThreeDisplay.Text  =  Format$(Total_Heat_rejected, 

"#,##0.0") 

History. CouplesStageThreeDisplay. Text  =  Format$(No_of_Couples, 

"#,##0.00") 

History.WperCoupleStageThreeDisplay.Text  =  Format$(P,  "#,##0.0000") 
History.TcStageThreeDisplay.Text  =  Format$(Tc,  "0.0000") 
History.ThStageThreeDisplay.Text  =  Format$(Th,  "0.0000") 
History.DeltaTStageThreeDisplay.Text  =  Format$(DeltaT,  "0.0000") 
History.  StageCOPStageThreeDisplay. Text  =  Format$(Coeff_of_Performance, 

"0.0000") 

History.OverallCOPStageThreeDisplay.Text  =  Format$(OverallCOP, 

"0.0000") 

History.Width  =  4425 
Case  4 

History.qLoadStageFourDisplay.Text  =  Format$(qT,  "#,##0.0") 
History. qRejectedStageFourDisplay.Text  =  Format$(Total_Heat_rejected, 

"#,##0.0") 

History. CouplesStageFourDisplay.Text  =  Format$(No_of_Couples, 

"#,##0.00") 

History. WperCoupleStageFourDisplay.Text  =  Format$(P,  "#,##0.0000") 
History.TcStageFourDisplay.Text  =  Format$(Tc,  "0.0000") 
History.ThStageFourDisplay.Text  =  Format$(Th,  "0.0000") 
History.DeltaTStageFourDisplay.Text  =  Format$(DeltaT,  "0.0000") 
History. StageCOPStageFourDisplay.Text  =  Format$(Coeff_of_Performance, 

"0.0000") 

History.OverallCOPStageFourDisplay.Text  =  Format$(OverallCOP,  "0.0000") 
History.Width  =  5520 
Case  5 

History.qLoadStageFiveDisplay.Text  =  Format$(qT,  "#,##0.0") 
History.qRejectedStageFiveDisplay.Text  =  Format$(Total_Heat_rejected, 

"#,##0.0") 

History. CouplesStageFiveDisplay. Text  =  Format$(No_of_Couples, 

"#,##0.00") 

History.WperCoupleStageFiveDisplay.Text  =  Format$(P,  "#,##0.0000") 
History.TcStageFiveDisplay.Text  =  Format$(Tc,  "0.0000") 
History.ThStageFiveDisplay.Text  =  Format$(Th,  "0.0000") 
History.DeltaTStageFiveDisplay.Text  =  Format$(DeltaT,  "0.0000") 
History.StageCOPStageFiveDisplay.Text  =  Format$(Coeff_of_Performance, 

"0.0000") 

History.OverallCOPStageFiveDisplay.Text  =  Format$(OverallCOP,  "0.0000") 
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History.  Width  =  6585 
End  Select 

History. Show  0 

ScrollUpdate  =  False 

LowTempText.Text  =  Format$(Th,  "##,##0.0") 

ScrollUpdate  =  True 

ReferLoadText.Text  =  Format$(Total_Heat_rejected,  "##0.0") 

LowTempUpdate  =  False 
Tempscroll. Value  =  (Tc  +  DeltaT)  *  10 
LowTempUpdate  =  True 

Command2. Caption  =  "Cascade" 
Screen.  MousePointer  =  0 

Cascading  =  False 
End  Sub 

Sub  TempScroll_Change  () 
If  ScrollUpdate  =  True  Then 
Tc  =  Val(LowTempText.Text) 
Th  =  Tempscroll. Value  /  10 
DeltaT  =  Th  -  Tc 

HighTempDisplay.Text  =  Format$(Th,  "##,#00.0") 
DeltaTDisplay.Text  =  Format$(DeltaT,  "##,#00.0") 
End  If 
End  Sub 

Sub  LowTempText_Change  () 
If  LowTempUpdate  =  True  Then 

Tc  =  Val(LowTempText.Text) 

DeltaTMax  =  .5  *  z  *  (Tc  A  2) 

TMax  =  DeltaTMax  +  Tc 

Tempscroll. Max  =  Tc  *  10  +  1 
Tempscroll. Min  =  TMax  *  10  -  1 
If  Stage  =  0  Then 

Tempscroll.  Value  =  Tempscroll. Min 
End  If 

Th  =  Tempscroll.  Value  / 10 
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If  Cascading  =  False  Then 

DeltaT  -  Th  -  Tc 
End  If 

TMinDisplay.Text  =  Format$(Tc,  "##,##0.0") 

TMaxDisplay.Text  =  Format$(TMax,  "##,##0.0") 

DeltaTDisplay.Text  =  Format$( DeltaT,  "##,##0.0") 
End  If 
End  Sub 

Sub  PerformanceCmd_Click  () 

Screen. MousePointer  =  1 1 

Main. ShowGraphName. Enabled  =  True 

PerformanceCmd.Caption  =  "Wait" 

Command2.  Visible  =  -1 

History. Hide 

Stage  =  0 

DeltaT  =  Th  -  Tc 

Th_plus_Tc  =  Th  +  Tc 

gamma  =  Sqr(l  +  (z  /  2)  *  (Th_plus_Tc)) 

Io  =  (alpha  *  DeltaT)  /  (R  *  (gamma  -  1)) 

Vo  =  (alpha  *  DeltaT)  +  Io  *  R 

P  =  Io  *  Vo 

qT  =  Val(ReferLoadText.Text) 

qP  =  alpha  *  Io  *  Tc 

qj  =  .5  *  (Io  A  2)  *  R 

qF  =  K  *  DeltaT 

qo  =  qP  -  (qJ  +  qF) 

No_of_Couples  =  qT  /  qo 

Total_Heat_rejected  =  qT  +  No_of_Couples  *  P 

Coeff_of_Performance  =  qo  /  P 

TMaxDisplay.Text  =  Format$(DeltaTMax  +  Tc,  "##,##0.0") 
TMinDisplay.Text  =  Format$(Tc,  "##,##0.0") 
DeltaTDisplay.Text  =  Format$(DeltaT,  "##,##0.0") 

Results.CoeffText.Text  =  Format$(Coeff_of_Performance,  "0.0000") 
Results.NumberCouplesText.Text  =  Format$(No_of_Couples,  "#,##0.00") 
Results.TotalHeatRejectedText.Text  =  Format$(Total_Heat_rejected,  "#,##0.0") 
Results.AperCoupleText.Text  =  Format$(Io,  "#,##0.00") 
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Results.VperCoupleText.Text  -  Format$(Vo,  "#,##0.0000") 

Results.WperCoupleText.Text  =  Format$(P,  "#,##0.00") 

Results. qTText.Text  =  Format$(qT,  "#,##0.0") 

Results.qPText.Text  =  Format$(qP,  "#,##0.0000") 

Results.qJText.Text  =  Format$(qJ,  "#,##0.0000") 

Results.qFText.Text  =  Format$(qF,  "#,##0.0000") 

Results.qoText.Text  =  Format$(qo,  "#,##0.0000") 

Results. Caption  =  "Optimum  Coefficient  of  Performance" 

Results.  Show 

Results.PowerConsumptionDisplay.Text  =  Format$((P  *  No_of_Couples), 
"#,###,##0.0") 

PerformanceCmd.  Width  =  975 

PerformanceCmd.Caption  =  "Restart" 

Screen. MousePointer  =  0 
End  Sub 

Cooler.txt 

//  Code  for  Optimum  Coeficient  of  Performance  Window 
DefDbl  A-Z 

Global  Const  True  =  -1 
Global  Const  False  =  0 
Global  Const  pi  =  355  /  113 
Global  DeltaT  As  Double 
Global  DeltaTMax  As  Double 

Global  alpha_A,  alpha_B,  rho_A,  rho_B,  kappa_A,  kappa_B  As  Double 

Global  alpha  As  Double 

Global  z,  K,  R  As  Double 

Global  Stage  As  Integer 

Global  ScrollUpdate  As  Integer 

Global  LowTempUpdate  As  Integer 

Global  Cascading  As  Integer 

Global  Total_Heat_rejected  As  Double 

Global  Tc  As  Double 

Global  Th  As  Double 

Global  TotalPower  As  Double 

Global  FirstStage_qo  As  Double 

Cascade.txt 

//  Code  for  Cascaded  Stages  Window 
Sub  PrintName_Click  () 
C  ascade .  PrintForm 
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End  Sub 

Sub  ExitName_Click  () 

Cascade.  Hide 
End  Sub 

Define.txt 

//  Code  to  call  Definitions  Window 
Sub  PrintName_Click  () 
Definitions.  PrintForm 
End  Sub 

Sub  ExitName_Click  () 

Definitions.Hide 
End  Sub 

Diagram.txt 

//  Code  to  call  Diagrams 
Sub  PrintName_Click  () 

Cooler.PrintForm 
End  Sub 

Sub  ExitName_Click  () 

Cooler.Hide 
End  Sub 

Ending.txt 

//  Code  to  end  design  session,  closing  all  open  windows 
Sub  Command  l_Click  () 

End 
End  Sub 

Sub  Command2_Click  () 

Ending. Hide 
End  Sub 

Sub  Form_Click  () 

Main.Left  =  (Screen.Width  -  Main.Width)  /  2 
Main. Top  =  (Screen. Height  -  Main. Height)  /  2 

End  Sub 

Sub  Form_Load  () 


56 


Ending.Left  =  (Screen.Width  -  Ending. Width)  /  2 
Ending. Top  =  (Screen. Height  -  Ending. Height)  /  2 

End  Sub 

Graph.txt 

//  Code  for  graping  routines 
Dim  Shown  As  Integer 
Dim  Samples  As  Integer 
Dim  IMax  As  Double 
Dim  IMin  As  Double 
Dim  YMax  As  Double 
Dim  YMin  As  Double 
Dim  OldYMax  As  Double 
Dim  OldlMax  As  Double 
Dim  OldlMin  As  Double 
Dim  OldMaxCount  As  Double 
Dim  OldMinCount  As  Double 
Dim  IatOldMinCount  As  Double 
Dim  IatOldMaxCount  As  Double 
Dim  PerformancePlot  As  Integer 
Dim  SuperimposePlots  As  Integer 

Sub  Graphlt  () 

ReDim  I(  Samples) 

ReDim  Y( Samples) 

IMin  =  ((2  *  alpha  *  Tc)  -  Sqr((2  *  alpha  *  Tc)  A  2  -  (8  *  R  *  K  *  DeltaT)))  /  (2  *  R) 

IMax  =  ((2  *  alpha  *  Tc)  +  Sqr((2  *  alpha  *  Tc)  A  2  -  (8  *  R  *  K  *  DeltaT)))  /  (2  *  R) 

YMin  =  0 

YMax  =  0 

MinCount  =  0 

MaxCount  =  0 

Minimum  =  False 

Steplnterval  =  (IMax  -  IMin)  /  (Samples  -  1) 

If  PerformancePlot  =  True  Then 

YAxisLabel.Text  =  "COP" 

Graph. Caption  =  "Performance" 
Else 

YAxisLabel.Text  =  "q" 

Graph. Caption  =  "Heat  Pumping" 
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End  If 

For  J  =  1  To  Samples 

I(J)  =  IMin  +  Steplnterval  *  (J  -  1) 

If  PerformancePlot  =  True  Then 

Y(J)  =  ((alpha  *  I(J)  *  Tc)  -  (.5  *  I(J)  A  2  *  R)  -  (K  *  DeltaT))  /  (I(J)  *  ((alpha  * 
DeltaT)  +  (I(J)  *  R))) 
Else 

Y(J)  =  (alpha  *  Tc  *  I(J))  -  (.5  *  I(J)  A  2  *  R)  -  (K  *  DeltaT) 
End  If 

If  Y(J)  >  0  And  Minimum  =  False  Then 

YMin  =  Y(J) 

MinCount  =  J 

Minimum  =  True 
ElseIfY(J)>YMaxThen 

YMax  =  Y(J) 

I_at_YMax  =  I(J) 
ElseIfY(J)>OThen 

MaxCount  =  J 
End  If 

Next  J 

If  SuperimposePlots  =  True  Then 

ResultsPicture.ScaleMode  =  0 

ResultsPicture.ScaleTop  =  Old  YMax 

ResultsPicture.ScaleHeight  =  -Old  YMax 

ResultsPicture.ScaleLeft  =  IatOldMinCount 

ResultsPicture.  Scale  Width  =  IatOldMaxCount  -  IatOldMinCount 

YMaxText.Text  =  Format$(01dYMax,  "0.00  e+00") 

IMinText.Text  =  Format$(01dIMin,  "0.00  e+00") 

IMaxText.Text  =  Format$(01dIMax,  "0.00  e+00") 
Else 

ResultsPicture. Cls 

ResultsPicture.ScaleMode  =  0 

ResultsPicture.ScaleTop  =  YMax 

ResultsPicture.ScaleHeight  =  -YMax 

ResultsPicture.ScaleLeft  =  I(MinCount) 

ResultsPicture. ScaleWidth  =  I(MaxCount)  -  I(MinCount) 

YMaxText.Text  =  Format$(YMax,  "0.00  e+00") 

IMinText.Text  =  Format$(IMin,  "0.00  e+00") 

IMaxText.Text  =  Format$(IMax,  "0.00  e+00") 
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OldYMax  =  YMax 
OldlMax  =  IMax 
OldlMin  =  IMin 
OldMaxCount  =  MaxCount 
IatOldMinCount  =  I(MinCount) 
IatOldMaxCount  =  I(MaxCount) 
End  If 

For  J  =  1  To  Samples 

ResultsPicture.Line  (I(J  -  1),  Y(J  -  1))-(I(J),  Y(J)) 

Next  J 

I_at_YMaxDisplay.Text  =  Format$(I_at_YMax,  "0.00  e+00") 
DeltaTDisplay.Text  =  Format$(DeltaT,  "0.00") 

End  Sub 

Sub  ExitName_Click  () 

Graph.  Hide 
End  Sub 

Sub  PrintName_Click  () 

Graph.  PrintForm 
End  Sub 

Sub  PerformanceName_Click  () 

PerformancePlot  =  True 

PerformanceName. Checked  =  True 

PumpName. Checked  =  False 
End  Sub 

Sub  PumpName_Click  () 

PerformancePlot  =  False 

PerformanceName. Checked  =  False 

PumpName. Checked  =  True 
End  Sub 

Sub  Form_Load  () 

PerformancePlot  =  True 

Samples  =  200 
End  Sub 
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Sub  DrawName_Click  () 

Screen. MousePointer  =  1 1 

Graphlt 

SuperName.Enabled  =  True 

Screen. MousePointer  =  0 
End  Sub 

SubalOOName_Click() 
a  lOOName.  Checked  =  True 
a200Name.  Checked  =  False 
a500Name.  Checked  =  False 
a  lOOOName.  Checked  =  False 
Samples  =  100 

End  Sub 

Sub  a200Name_Click  () 
a  lOOName.  Checked  =  False 
a200Name.  Checked  =  True 
a500Name. Checked  =  False 
a  lOOOName. Checked  =  False 
Samples  =  200 

End  Sub 

Sub  a500Name_Click  () 
a  lOOName.  Checked  =  False 
a200Name.Checked  =  False 
a500Name.  Checked  =  True 
a  lOOOName. Checked  =  False 
Samples  =  500 

End  Sub 

SubalOOOName_Click() 
alOOName. Checked  =  False 
a200Name.Checked  =  False 
a500Name. Checked  =  False 
a  lOOOName. Checked  =  True 
Samples  =  1000 

End  Sub 

Sub  SuperName_Click  () 

If  SuperName. Checked  =  False  Then 
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SuperimposePlots  =  True 
SuperName. Checked  =  True 

Else 

SuperimposePlots  =  False 
SuperName. Checked  =  False 

End  If 
End  Sub 

History.txt 

//  Cascaded  Stages  Window  menu  bar 
Sub  Form_Load  () 

History.Left  =  0 

History. Top  =  Main. Height 
End  Sub 

Sub  PrintName_Click  () 

History.PrintForm 
End  Sub 

Intro.txt 

//  Introduction  graphic,  program  startup 
Sub  Timer  l_Timer  () 

Screen. MousePointer  =  1 1 

Load  Results 

Load  Performance 

Load  Material 

Load  Definitions 

Load  TempCalc 

Load  Graph 

Load  Cascade 

Main. Show 

Load  Cooler 

Load  Ending 

Load  History 

Intro. Hide 

Timerl.  Enabled  =  0 

Screen. MousePointer  =  0 
End  Sub 

Sub  Form_Load  () 

Screen. MousePointer  =  1 1 
End  Sub 
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Pumping.txt 

//  Code  for  Maximum  Heat  Pumped  Window 
Sub  Command l_Click  () 

Screen.  MousePointer  =  1 1 

Stage  =  0 

Pumping. Hide 

History. Hide 

Results. Hide 

Screen. MousePointer  =  0 
End  Sub 

Sub  Form_Load  () 

PerformanceCmd.  Width  =  2415 

Pumping. Left  =  Screen. Width  -  Pumping. Width 

Pumping. Top  =  Main. Height  +  450 

LowTempUpdate  =  True 

Cascading  =  False 

End  Sub 

Sub  Command2_Click  () 

Screen.  MousePointer  =  1 1 
Command2.  Caption  =  "Wait" 
History.WindowState  =  0 
LowTempUpdate  =  True 
Cascading  =  True 

Stage  =  Stage  +  1 

If  Stage  =  1  Then 
TotalPower  =  0 
End  If 

Tc  =  Val(LowTempText.Text) 
Th  =  Val(HighTempDisplay.Text) 

DeltaT  =  Th  -  Tc 

DeltaTMax  =  .5  *  z  *  (Tc  A  2) 

Th_plus_Tc  =  Th  +  Tc 

gamma  =  Sqr(l  +  (z  /  2)  *  (Th_plus_Tc)) 

Io  =  (alpha  *  Tc)  /  R 
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Vo  =  alpha  *  Th 

P  =  Io  *  Vo 

qT  =  Val(ReferLoadText.Text) 

qP  =  alpha  *  Io  *  Tc 

qj  =  .5  *  (Io  A  2)  *  R 

qF  =  K  *  DeltaT 

qo  =  .5  *  ((alpha  A  2  *  Tc  A  2)  /  R)  -  qF 

No_of_Couples  =  qT  /  qo 

Total_Heat_rejected  =  qT  +  No_of_Couples  *  P 

Coeff_of_Performance  =  qo  /  P 

TotalPower  =  TotalPower  +  (P  *  No_of_Couples) 

If  Stage  =  1  Then 

FirstStage_qo  =  qo  *  No_of_Couples 
End  If 
OverallCOP  =  FirstStage_qo  /  TotalPower 

TMaxDisplay.Text  =  Format$(DeltaTMax  +  Tc,  "##,##0.0") 
TMinDisplay.Text  =  Format$(Tc,  "##,##0.0") 
DeltaTDisplay.Text  =  Format$(DeltaT,  "##,##0.0") 

Results.CoeffText.Text  =  Format$(Coeff_of_Performance,  "0.0000") 
Results.NumberCouplesText.Text  =  Format$(No_of_Couples,  "#,##0.00") 
Results.TotalHeatRejectedText.Text  =  Format$(Total_Heat_rejected,  "#,##0.0") 
Results.AperCoupleText.Text  =  Format$(Io,  "#,##0.00") 
Results.VperCoupleText.Text  =  Format$(Vo,  "#,##0.0000") 
Results.WperCoupleText.Text  =  Format$(P,  "#,##0.0000") 
Results.qTText.Text  =  Format$(qT,  "#,##0.0") 
Results.qPText.Text  =  Format$(qP,  "#,##0.0000") 
Results.qJText.Text  =  Format$(qJ,  "#,##0.0000") 
Results.qFText.Text  -  Format$(qF,  "#,##0.0000") 
Results.qoText.Text  =  Format$(qo,  "#,##0.0000") 
Results.Caption  =  "Maximum  Heat  Pumping"  +  "  Stage  "  +  Str$(Stage) 
Results.  Show 

Select  Case  Stage 
Case  1 

History.qLoadStageOneDisplay.Text  =  Format$(qT,  "#,##0.0") 
History.qRejectedStageOneDisplay.Text  =  Format$(Total_Heat_rejected, 

"#,##0.0") 

History. CouplesStageOneDisplay.Text  =  Format$(No_of_Couples, 

"#,##0.00") 

History. WperCoupleStageOneDisplay .Text  =  Format$(P,  "#,##0.0000") 
History.TcStageOneDisplay.Text  =  Format$(Tc,  "0.0000") 
History.ThStageOneDisplay.Text  =  Format$(Th,  "0.0000") 
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History.DeltaTStageOneDisplay.Text  =  Format$(DeltaT,  "0.0000") 
History. StageCOPStageOneDisplay. Text  =  Format$(Coeff_of_Performance, 

"0.0000") 

History.OverallCOPStageOneDisplay.Text  =  Format$(OverallCOP,  "0.0000") 
History.Width  =  2250 
Case  2 

History.qLoadStageTwoDisplay.Text  =  Format$(qT,  "#,##0.0") 
History.  qRejectedStageTwoDisplay.Text  =  Format$(Total_Heat_rejected, 

"#,##0.0") 

History. CouplesStageTwoDisplay. Text  =  Format$(No_of_Couples, 

"#,##0.00") 

History.WperCoupleStageTwoDisplay.Text  =  Format$(P,  "#,##0.0000") 
History.TcStageTwoDisplay.Text  =  Format$(Tc,  "0.0000") 
History.ThStageTwoDisplay.Text  =  Format$(Th,  "0.0000") 
History.DeltaTStageTwoDisplay.Text  =  Format$(DeltaT,  "0.0000") 
History.  StageCOPStageTwoDisplay. Text  =  Format$(Coeff_of_Performance, 

"0.0000") 

History.OverallCOPStageTwoDisplay.Text  =  Format$(OverallCOP,  "0.0000") 
History.Width  =  3360 
Case  3 

History.qLoadStageThreeDisplay.Text  =  Format$(qT,  "#,##0.0") 
History.  qRejectedStageThreeDisplay.Text  =  Format$(Total_Heat_rejected, 

"#,##0.0") 

History. CouplesStageThreeDisplay.Text  =  Format$(No_of_Couples, 

"#,##0.00") 

History. WperCoupleStageThreeDisplay.Text  =  Format$(P,  "#,##0.0000") 
History.TcStageThreeDisplay.Text  =  Format$(Tc,  "0.0000") 
History.ThStageThreeDisplay.Text  =  Format$(Th,  "0.0000") 
History.DeltaTStageThreeDisplay.Text  =  Format$(DeltaT,  "0.0000") 
History.  StageCOPStageThreeDisplay. Text  =  Format$(Coeff_of_Performance, 

"0.0000") 

History.OverallCOPStageThreeDisplay.Text  =  Format$(OverallCOP, 

"0.0000") 

History.Width  =  4425 
Case  4 

History.qLoadStageFourDisplay.Text  =  Format$(qT,  "#,##0.0") 
History.qRejectedStageFourDisplay.Text  =  Format$(Total_Heat_rejected, 

"#,##0.0") 

History. CouplesStageFourDisplay.Text  =  Format$(No_of_Couples, 

"#,##0.00") 

History.WperCoupleStageFourDisplay.Text  =  Format$(P,  "#,##0.0000") 
History.TcStageFourDisplay.Text  =  Format$(Tc,  "0.0000") 
History.ThStageFourDisplay.Text  =  Format$(Th,  "0.0000") 
History.DeltaTStageFourDisplay.Text  =  Format$(DeltaT,  "0.0000") 
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History.  StageCOPStageFourDisplay.  Text  =  Format$(Coeff_of_Performance, 

"0.0000") 

History.OverallCOPStageFourDisplay.Text  =  Format$(OverallCOP,  "0.0000") 
History.Width  -  5520 
Case  5 

History.qLoadStageFiveDisplay.Text  =  Format$(qT,  "#,##0.0") 
History. qRejectedStageFiveDisplay. Text  =  Format$(Total_Heat_rejected, 

"#,##0.0") 

History.CouplesStageFiveDisplay.Text  =  Format$(No_of_Couples, 

"#,##0.00") 

History.WperCoupleStageFiveDisplay.Text  =  Format$(P,  "#,##0.0000") 
History.TcStageFiveDisplay.Text  =  Format$(Tc,  "0.0000") 
History.ThStageFiveDisplay.Text  =  Format$(Th,  "0.0000") 
History.DeltaTStageFiveDisplay.Text  =  Format$(DeltaT,  "0.0000") 
History. StageCOPStageFiveDisplay. Text  =  Format$(Coeff_of_Performance, 

"0.0000") 

History.OverallCOPStageFiveDisplay.Text  =  Format$(OverallCOP,  "0.0000") 
History.Width  =  6585 
End  Select 

History. Show  0 

ScrollUpdate  =  False 

LowTempText.Text  =  Format$(Th,  "##,##0.0") 

ScrollUpdate  =  True 

ReferLoadText.Text  =  Format$(Total_Heat_rejected,  "##0.00") 

LowTempUpdate  =  False 
Tempscroll. Value  =  (Tc  +  DeltaT)  *  10 
LowTempUpdate  =  True 

Command2. Caption  =  "Cascade" 
Screen.  MousePointer  =  0 

Cascading  =  False 
End  Sub 

Sub  TempScroll_Change  () 
If  ScrollUpdate  =  True  Then 
Tc  =  Val(LowTempText.Text) 
Th  =  Tempscroll.  Value  / 10 
DeltaT  =  Th  -  Tc 
HighTempDisplay.Text  =  Format$(Th,  "##,#00.0") 
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DeltaTDisplay.Text  =  Format$(DeltaT,  "##,#00.0") 
End  If 
End  Sub 

Sub  LowTempText_Change  () 
If  LowTempUpdate  =  True  Then 

Tc  =  Val(LowTempText.Text) 

DeltaTMax  =  .5  *  z  *  (Tc  A  2) 

TMax  =  DeltaTMax  +  Tc 

Tempscroll.Max  =  Tc  *  10  +  1 
Tempscroll.Min  =  TMax  *  10  -  1 
If  Stage  =  0  Then 

Tempscroll.  Value  =  Tempscroll.Min 
End  If 

Th  =  Tempscroll. Value  /  10 

If  Cascading  =  False  Then 

DeltaT  =  Th  -  Tc 
End  If 

TMinDisplay.Text  =  Format$(Tc,  "##,##0.0") 
TMaxDisplay.Text  =  Format$(TMax,  "##,##0.0") 
DeltaTDisplay.Text  =  Format$( DeltaT,  "##,##0.0") 

End  If 

End  Sub 

Sub  PerformanceCmd_Click  () 

Screen.  MousePointer  =  1 1 

Main.  ShowGraphName.  Enabled  =  True 

PerformanceCmd.  Caption  =  "Wait" 

Command2.  Visible  =  -1 

History. Hide 

Stage  -  0 

DeltaT  =  Th  -  Tc 

Th_plus_Tc  =  Th  +  Tc 

gamma  =  Sqr(l  +  (z  /  2)  *  (Th_plus_Tc)) 

Io  =  (alpha  *  Tc)  /  R 
Vo  =  alpha  *  Th 
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P  =  Io  *  Vo 

qT  =  Val(ReferLoadText.Text) 

qP  =  alpha  *  Io  *  Tc 

qj  =  .5  *  (Io  A  2)  *  R 

qF  =  K  *  DeltaT 

qo  =  .5  *  ((alpha  A  2  *  Tc  A  2)  /  R)  -  qF 

No_of_Couples  =  qT  /  qo 

Total_Heat_rejected  =  qT  +  No_of_Couples  *  P 

Coeff_of_Performance  =  qo  /  P 


TMaxDisplay.Text  =  Format$(DeltaTMax  +  Tc,  "##,##0.0") 
TMinDisplay.Text  =  Format$(Tc,  "##,##0.0") 
DeltaTDisplay.Text  =  Format$(DeltaT,  "##,##0.0") 

Results.CoeffText.Text  =  Format$(Coeff_of_Performance,  "0.0000") 

Results.NumberCouplesText.Text  =  Format$(No_of_Couples,  "#,##0.00") 

Results.TotalHeatRejectedText.Text  =  Format$(Total_Heat_rejected,  "#,##0.0") 

Results.AperCoupleText.Text  =  Format$(Io,  "#,##0.00") 

Results. VperCoupleText.Text  =  Format$(Vo,  "#,##0.0000") 

Results. WperCoupleText.Text  =  Format$(P,  "#,##0.00") 

Results.qTText.Text  =  Format$(qT,  "#,##0.0") 

Results.qPText.Text  =  Format$(qP,  "#,##0.0000") 

Results.qJText.Text  =  Format$(qJ,  "#,##0.0000") 

Results.qFText.Text  =  Format$(qF,  "#,##0.0000") 

Results.qoText.Text  =  Format$(qo,  "#,##0.0000") 

Results. Caption  =  "Maximum  Heat  Pumping" 

Results.  Show 

Results. PowerConsumptionDisplay.Text  =  Format$((P  *  No_of_Couples), 
"#,###,##0.00") 

PerformanceCmd.Width  =  975 

PerformanceCmd.Caption  =  "Restart" 

Screen.  MousePointer  =  0 
End  Sub 


Results.txt 

//  Code  for  Results  Window 

Sub  Form_Click  () 

Results.Top  =  (Screen. Height  -  Results. Height) 
Results.Left  =  (Screen.Width  -  Results. Width)  /  2 

End  Sub 
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Sub  Form_Load  () 

Results. Top  =  (Screen. Height  -  Results.Height) 
Results.Left  =  (Screen. Width  -  Results. Width)  /  2 

End  Sub 

Sub  PrintName_Click  () 

Results. PrintForm 
End  Sub 

Tempcalc.txt 

//  Code  for  Temperature  Conversion  Calculator 
Sub  Command2_Click  () 

ConvertDisplayText.Text  =  "" 

ConvertDisplayText.SetFocus 
End  Sub 

Sub  ConvertCmd_Click  () 
Const  True  =  -1 
Const  False  =  0 

Temperature  =  Val(ConvertDisplayText.Text) 
If  FarenheitOutOption.  Value  =  True  Then 
If  RankinelnOption.  Value  =  True  Then 

Temperature  =  Temperature  -  459.67 
Elself  FarenheitlnOption.  Value  =  True  Then 

Temperature  =  Temperature 
Elself  CelsiusInOption.  Value  =  True  Then 

Temperature  =  (Temperature  *  9  /  5)  +  32 
Elself  KelvinlnOption. Value  =  True  Then 

Temperature  =  ((Temperature  -  273.15)  *  9  /  5)  +  32 
End  If 

Elself  RankineOutOption.  Value  =  True  Then 

If  RankinelnOption.  Value  =  True  Then 

Temperature  =  Temperature 
Elself  FarenheitlnOption.  Value  =  True  Then 

Temperature  =  Temperature  +  459.67 
Elself  CelsiusInOption.  Value  =  True  Then 

Temperature  =  ((Temperature)  *  9  /  5)  +  32  +  459.67 
Elself  KelvinlnOption.  Value  =  True  Then 

Temperature  =  ((Temperature  -  273.15)  *  9  /  5)  +  32  +  459.67 
End  If 
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Elself  CelsiusOutOption. Value  =  True  Then 

If  RankinelnOption.  Value  =  True  Then 

Temperature  =  (Temperature  -  32  -  459.67)  *  5  /  9 
Elself  FarenheitlnOption. Value  =  True  Then 

Temperature  =  (Temperature  -  32)  *  5  /  9 
Elself  CelsiusInOption.  Value  =  True  Then 

Temperature  =  Temperature 
Elself  KelvinlnOption. Value  =  True  Then 

Temperature  =  Temperature  -  273.15 
End  If 

Elself  KelvinOutOption. Value  =  True  Then 

If  RankinelnOption.  Value  =  True  Then 

Temperature  =  (Temperature  -  32  -  459.67)  *  5  /  9  +  273.15 
Elself  FarenheitlnOption.  Value  =  True  Then 

Temperature  =  (Temperature  -  32)  *  5  /  9  +  273.15 
Elself  CelsiusInOption.  Value  =  True  Then 

Temperature  =  Temperature  +  273.15 
Elself  KelvinlnOption.  Value  =  True  Then 

Temperature  =  Temperature 
End  If 


End  If 

ConvertDisplayText.Text  =  Format$(Temperature,  "##,##0.00") 

ConvertDisplayText.SetFocus 

End  Sub 

Sub  Form_Click  () 

ConvertDisplayText.SetFocus 
End  Sub 

Sub  Form_Load  () 

TempCalc.Left  =  (Screen.Width  -  TempCalc.Width)  /  2 
TempCalc.Top  =  (Screen. Height  -  TempCalc. Height)  /  2 

End  Sub 

Sub  ExitName_Click  () 

TempCalc. Hide 
End  Sub 
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